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ABtTlMCT 


tosoMnt  fraquanoy  aqiwtlona  for  right  circular  cylindora  of  any 
diBonaiona  wora  darivad  to  aid  in  tha  daaign  of  Bicrominiatura  caraaic' 
piaBoalactrio  filtara  and  tranafonaara .  Tha  claaaical  Raylaigh  nathod 
vaa  BOdifiad  to  includa  a  phantoa  work  tana  that  accountad  for  anar- 
giaa  loat  to  vibrationa  in  tha  eouplad  Boda.  Tha  raaulting  aquationa 
wara 

...... .(i;  ^ 

^  Lomtudiui  »d.  (sty  ♦  (|7  I II  ]  -  I 


whara  f  and  a  ara  tha  fundaaantal  fraquanciaa;  0  and  H  ara  fraquancy 

conatanta;  L  and  D  ara  langth  ami  dianatar;  and  ara  Poiaaon'a 

ration;  and  R  in  a  Baaaal  function  ralatad  to  Poiaaonia  ratio, 
o 

Tha  aquationa  darivad  wara  varifiad  axpariBantally  ovar  a  full 
ranga  of  dlBonaion  ratloa  for  BaTiOg  cylindara  and  for  a  nnall  group 
of  ataal  cyllndam. 

1.  IKTHODUCnOW 


Kaact  aquationa  for  datamining  tha  raaonant  fraquanclaa  of  cy- 
llndrically  abapad  bodlaa  ara  avallabla  in  the  literature  for  ideal 
one-diaenalonal  caaaa,  i.e.  for  vary  thin  diaka  and  anall  diawatar 
roda  (raf  1). 

'o  ■ 

(2) 

f  ■  Idaal  radial  raaonant  fraquancy  for  diak  with  aaro  thicknaaa,  kc 
.  o 

d  ■  Idaal  longitudinal  raaonant  fraquancy  for  rod  with  aero  dlaBOter,  kc 
o 

D  ■  DiaBoter,  In. 

L  •  Langth,  in. 

I  m  Radial  Bodulua  of  alaaticlty,  Ib/ln.* 

R 


Long  Rod 


BLV  p 


S 


E-  3  Longitudinal  modulus  of  elasticity,  Ib/ln.^ 

la 

p  3  Density,  lb“sec*/in.* 

3  Poisson's  ratio  associating  radial  and  tangential  strains  in  the 
radial  plane,  dimensionless 

J  3  Bessel  function 


The  tiny  size  of  microminiature  piezoelectric  elements  usually  makes 
the  second  dimension  significant  and  the  ideal  handbook  equations  then 
become  inapplicable.  Attempts  reported  in  the  literature  to  take  into 
account  the  effects  of  the  second  dimension  are  adequate  only  over  a 
narrow  range  of  lengtb-diameter  ratios. 

Lord  Rayleigh  (ref  2)  applied  the  principle  of  conservation  of  ener¬ 
gy  to  a  long,  thin  isotropic  cylinder  to  derive  the  resonant  frequency 
equation  for  the  longitudinal  mode.  He  accounted  for  the  kinetic  energy 
in  the  radial  mode  during  longitudinal  resonance  by  assuming  that 
Poisson's  ratio  may  be  used  to  describe  the  radial  deflections  in  terms 
of  the  longitudinal.  He  noted  that  the  inertia  of  the  particles  was  dis¬ 
regarded  in  making  the  assumption,  but  felt  it  would  be  adequate  for  very 
small  diameters.  Inclusion  of  the  radial  energy  in  the  equation  for  con¬ 
servation  of  energy  modified  the  ideal  longitudinal  frequency  equation 
to: 


Mosley  (ref  3)  used  Newton's  laws  and  the  theory  of  elasticity  to 
determine  the  effect  of  length  on  the  radial  vibrations  of  a  thick  disk. 
However,  to  obtain  a  solution,  simplifying  assumptions  were  necessary 
and  limited  the  use  of  his  equation 


to  small  length-diameter  ratios. 

In  this  report,  the  longitudinal  and  radial  frequency  equations  are 
derived  for  cylinders  of  any  dimensions.  Experimental  results  from 
piezoelectric  ceramic  cylinders  are  then  used  to  verify  the  equations. 

2.  THEORY 


The  modified  Rayleigh  procedure  Incorporates  the  following  steps: 

(1)  Imagine  that  energy  for  radial  vibrations  is  supplied  by  de¬ 
forming  the  cylinder  slowly  with  an  external  pressure  acting  radially. 
Calculate  the  work  done  by  this  pressure. 
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(2)  Find  the  energy  lost  to  the  thickness  (longitudinal)  mode 
by  defining  a  "phantom”  work  as  that  work  which  would  have  resulted 
had  the  only  existing  force  been  a  phantom  force  compressing  the 
cylinder  length.  The  phantom  force  as  Interpreted  in  this  definition 
is  the  product  of  some  constant  N|^  and  the  net  radial  force. 

(3)  Specify  the  work  available  for  radial  vibrations  by  subtract¬ 
ing  the  phantom  work  from  the  radial  work.  Assuming  that  the  radial 
motion  of  the  particles  will  be  sinusoidal  when  the  pressure  is  suddenly 
released,  an  expression  for  the  radial  resonant  frequency  as  a  function 
of  dimensions  may  be  derived  by  equating  the  available  radial  work  to 
the  maximum  uncoupled  kinetic  energy. 

(4)  An  expression  for  longitudinal  resonant  frequency  as  a  func¬ 
tion  of  dimensions  is  derived  by  equating  the  maximum  kinetic  energy 
in  the  lengthwise  plane  and  the  difference  between  the  work  initially 
supplied  to  the  length  and  a  phantom  work  acting  radially.  The  phan¬ 
tom  force  in  this  case  is  the  product  of  the  net  longitudinal  force 
and  some  constant  N,  . 

li 

The  conditions  imposed  upon  the  derivation  to  follow  are: 

Right  circular  cylinder 
Free,  undamped  vibration 
Fundamental  frc^quency 
Homogeneous  material 
Planar  isotropy 
Linearly  elastic  material 
Constant  temperature 

The  resonant  frequency  equations  of  a  right  circular  cylinder  of 
any  length  and  diameter  will  now  be  derived  using  a  simple  spring-mass 
analogy  and  the  procedure  outlined  above.  The  equivalent  spring  con¬ 
stants  and  ideal  frequencies  in  the  derivation  are  presented  below  for 
the  extreme  cases  of  the  cylinder,  namely  the  thin  disk  and  the  long 
rod. 


2.1  Spring-Mass  Analogy  of  a  Long,  Thin  Rod  for  Free  Lengthwise 
Vibration 


F. 


I  L  I 


'  6.  I 


vm- 


/ 

/  K 

A  h 


Fj^(L/2) 

®L  =  FW/l) 
F  E  Trt)S 

K  -  —  - 
^  ®L 


2L 


a  = 
o 


^1/"l 


n  1  . 1  jj  n 
u  =  2L”  *  “  * 


(5) 


(6) 


(ref  4,  p  42) 
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6,  a  Maximum  longitudinal  deflection.  In. 

F,  a  Net  longitudinal  force,  lb 
L 

K  a  Longitudinal  spring  constant,  Ib/ln. 

li 

a  Effective  lumped  mass  of  cylinder  In  longitudinal  vibration, 
lb-8ec*/in. 

H  a  Longitudinal  frequency  constant,  Inrrkc 

2.2  Spring-Mass  Analogy  of  a  Thin  Disk  for  Free  Radial  Vibration 


I 


Zn  h 


V  !r  ^ 
»  °  ‘r  ■ 


1  • 
16-  ‘ 
I  R  ' 


j./!" 

2Tr\/  Mj^  “  W)]  p(l-p 


(Ref  4,  p  306) 


6-  a  Maximum  radial  deflection.  In.  ’ 

R 

K-  a  Radial  spring  constant,  Ib/ln, 

R 

Fj^  a  Net  radial  force,  lb 
0  a  Radial  frequency  constant,  in.-kc 

M_  a  Effective  lumped  mass  of  cylinder  in  radial  vibration,  Ib-sec^/ln, 
R 

Effective  values  of  E  and  p  are  used  for  anisotropic  materials. 

For  plesoelectrlcs,  the  effective  values  are  those  that  would  be  ob¬ 
tained  experimentally  with  the  sample  fully  polarised. 

3.  DERIVATION  OF  BOIATIONS 


Now  that  the  equivalent  spring  systems  and  their  properties  are 
defined,  the  generalized  frequency  equations  may  be  derived  In  terms  of 
these  systems.  The  radial  mode  will  be  attacked  first 

(1)  Deteivlna  the  work  done 'bp  the  cylinder  b|f  a  tensile  pressure 
acting  in  the  radial  direction  .  ; 


w  •  ^  p  A  ^  R 

*ln  2  ^R  ®R  2 
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(2)  Calculate  the  phantom  work  lost  to  vibrations  along  the 
length.  The  phantom  force  was  defined  as 


P  HR 


•  (12) 


Radial  constant  of  the  phantom  force 

Thus  the  phantom  work  is  given  by 


2  *P  “  2 


1 

2  K. 


(13) 


(3)  find  the  resonant  frequency  for  radial  vibrations 

(a)  Find  the  work  available  for  radial  vibrations 


2  it. 


(»,  F,)« 


(14) 


(b)  Find  the  maximum  uncoupled  kinetic  energy.  This  is 


found  from 


KK  »  ^  MV® 
max  2  max 


(15) 


V  =  Maximum  velocity  of  M„ 
max  '  R 

where  the  maximum  velocity  is  found  by  differentiating  the  uncoupled 
spring  displacement  =  6^^  cos  (Zirft)  with  respect  to  time  and  making 
the  result  a  maximum. 


When  this  is  done  the  maximum  kinetic  energy  becomes 

s 


KE 


max 


=  i  ■<^)  <*«>■ 


(16) 


(c)  Equate  the  energies  to  find  the  frequency  expression 


AW  a  KE 


max 


(17) 


(18) 


Introducing  the  spring  properties  (6,  9,  10)  equation 
(18)  simplifies  to 
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(19) 


where 


"U/ 


(20) 


Cl^ 


V  ®*R 


(21) 


(4)  An  analysis  slnllar  to  that  performed  for  the  radial  mode 
will  show  that  the  longitudinal  frequency  equation  Is  of  the  form 


(22) 


N  a  Longitudinal  constant  of  the  phantom  force. 
L 


becomes 


where 


Introducing  the  spring  properties ^6 ,  1,  9),  the  equation 


- 


^L 

4o 

V  »L 
‘  8 


(23) 

(24) 

(25) 


It  Is  apparent  that  the  radial  and  longitudinal  equations  (19,23) 
are  Incomplete  mathematically  In  their  present  form  In  that  they  con¬ 
tain  undefined  constants,  N-  and  N  .  If  left  In  this  form  the  equa¬ 
tions  would  be  purely  experimental^  thus,  It  Is  necessary  to  break  the 
constants  Into  the  more  familiar  material  properties.  This  may  be  ac¬ 
complished  by  observing  the  behavior  of  coupled  modes  Just  before  free 
vibration  occurs,  and  then  comparing  the  behavior  with  that  during 
free  vibration. 


A  radial  force  develops  a  deflection  in  the  longitudinal 
direction  due  to  Poisson  coupling.  The  magnitude  of  the  deflection 
of  each  end  of  the  cylinder  Is 


8 


L 


!ki: 

3 


^RL 

TT  D  K„ 


(26) 
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e.  m  Longitudinal  strain 
L 

Poisson's  ratio  relating  lengthwise  and  radial  strains. 

A  longitudinal  force  that  would  produce  the  same  deflection  Is: 


F 


L 


'l  'l 
2L 


(6) 


Consequently  the  static  relationship  between  the  lengthwise  and  the 
radial  force  Is: 


F 


L 


^  »  “l 

7^1 


(27) 


whereas  It  was  assumed  at  the  beginning  of  the  derivation  that  the  dy¬ 
namic  relationship  was: 


R  R 


(12) 


As  a  result j  the  phantom  force  will  be  expressed  as: 


f  ^ 


^R*R 


B  numerical  constant  (gg) 


Thus^  using  equation  (21) 


(29) 


and 


*R 


(30) 


when  the  arbitrary  assumption  that  Y?  ■  2  Is  made.  Verification  will 
be  done  experimentally. 
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Ordinarily  the  quantities  Bu  and  B|^  are  unavailable  In  the  literature. 
However,  their  ratio  is  proportional  to  the  frequency  constants  0  and 
H,  which  are  easily  obtained  from  literature  or  frequency  measurements 
on  a  long  rod  and  thin  disk.  should  therefore  be  expressed  In  terms 
of  the  frequency  constants  by  using  equations  20  and  24. 

4u  R 
0  IT 

A  similar  analysis  will  show  the  longitudinal  constant  to  be 


T  “  2  H  R 

o 


(33) 


Defining  these  constants  completely  generalises  the  frequency 
equations  19  and  23  to  functions  of  dimensions  and  material 

(34) 


(36) 

4.  VERIFICATION  OF  BqUATlOWB 
4.1  General  Considerations 


C^T<S)‘ 

"  ^RL  «  “o  1 

s 

«  1 

( 

[fhilJ 

\i.  0  ’T 

e 

-  1 

2  H  R 

•  o  ^ 

The  equations  were  verified  experlmentaliy  by  using  pleso- 
electrlc  ceramic  cylinders.  Being  plesoelectric  (ref  7,  p.  80),  means 
the  electrons  and  nuclei  within  t)w  ceramic  are  directionally  forced 
by  the  act  of  polarisation  Into  geometrical  positions  that  contain  no 
center  of  symmetry.  Consequently,  the  ceramic  assumes  a  motor  or  gen¬ 
erator  capability  In  that  applying  an  electric  field  causes  the  body 
to  expand  or  contract  according  to  the  direction  of  the  field,  and  con¬ 
versely,  physical  distortions  cause  the  body  to  develop  an  electric 
charge . 


Utilising  the  motor  ability,  the  body  may  be  excited  Into 
mechanical  vibration  by  an  alternating  electric  field.  If  the  field 
Is  applied  at  the  natural  resonant  frequency  of  the  body,  minimum  en¬ 
ergy  Is  required  to  Induce  maximum  amplitude  of  vibration.  The  gener¬ 
ator  action  of  the  plesoelectric  being  sensitive  to  the  amplitude 
of  strain  Is  electrically  a  minimum  Impedance  at  resonance.  As  a 
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result,  by  noting  which  frequency  Induces  s  minimum  Impedance,  the 
natural  resonant  frequency  of  the  body  Is  determined. 

Although  the  piezoelectric  cylinders  are  forced  by  an  al¬ 
ternating  electric  field,  the  natural  resonance  Is  Identical  with 
that  of  the  body  In  free  vibration.  This  Is  true  only  as  long  as 
there  are  no  external  mechanical  or  electrical  loads  on  the  device. 


Ilie  experimental  procedure  was  to  form  a  long  piezoelectric 
ceramic  cylinder,  and  then  record  the  resonance  each  time  part  of  the 
length  was  ground  away.  However,  nothing  other  than  the  strongest  funda- 
Biental  was  recorded  because  of  the  confusion  In  distinguishing  the  funda¬ 
mental  of  the  weak  coupled  mode  from  harmonics  and  subharmonics. 

Hm  experimental  results  were  used  to  show  that  the  form  of  the 
equatloni)ls  correct.  Hie  material  constants  were  then  obtained  experi¬ 
mentally  from  one  sample  and  used  to  show  that  they  are  valid  for  other 
samples  of  different  dimensions  but  of  the  same  material.  Finally,  data 
were  extracted  from  published  literature  and  compared  directly  with  the 
form  of  the  theoretical  equations. 

4.2  Bxperlmental  Methods,  Data,  and  Graphs 

Three  cylinders  were  made  of  a  mixture  of  barium  tltanate  and 
4  percent  lead  tltanate  and  were  cut  with  an  ultrasonic  Impact  grinder 
In  the  following  Initial  dimensions. 


Sample 

Diameter 

Length 

(In.) 

(in.) 

BT-1 

0.137 

0.307 

BT-2 

0.248 

0.159 

BT-3 

0.147 

0.359 

The  ends  were  silvered  with  alr-drylng  silver  paint,  and  then  the  sam¬ 
ples  were  polarized  longitudinally  by  Immersing  In  silicone  oil  heated 
to  120°C,  applying  40  v/mll  across  the  samples,  and  allowing  the  oil  to 
cool  with  voltage  maintained.  When  the  temperature  of  the  oil  bath 
dropped  to  about  50°C,  the  voltage  was  turned  off,  and  the  specimens  were 
removed,  washed,  and  allowed  to  cool  to  room  temperature. 

The  samples  were  placed  In  a  holder  consisting  of  two  leaf 
springs  having  a  very  low  mechanical  resonant  frequency.  Electrical  point 
coptact  was  made  to  the  silvered  surfaces  of  the  sample  under  light  pres¬ 
sure  near  the  center  of  the  disk.  Resonant  frequency  determinations  were 
made  at  room  temperature  with  a  Hewlett-Packard  Type  60eA  signal  generator 
and  a  Boonton  Electronics  lype  91B  RF  Voltmeter  as  Illustrated  In  figure  1. 
All  resonant  frequency  readings  were  made  Instantly  to  avoid  the  problems 
of  possible  Internal  heating  and  frequency  drift. 


13 


PIEZOELECTRIC  UNDER 

test 


i 

i 

Figure  1(a)  Apparatus  employed  for  measuring  the 
resonant  frequencies  of  piezoelectric  cylinders. 


Figure  l(h)  Leaf  spring  sample  holder. 
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The  lowest  frequency  producing  a  major  pip  on  the  voltmeter 
was  recorded  for  each  sample.  The  major  pip  Indicates  an  Impedance 
change  that,  when  related  to  the  motor  and  generator  action  of  the 
piezoelectric,  depicts  the  fundamental  mechanical  resonance. 

The  samples  were  then  ground  to  new  lengths,  the  ends  re- 
silvered,  and  the  frequency  measurements  repeated.  The  resonances  ob¬ 
tained  are  shown  In  table  I . 

Table  II  lists  the  data  extracted  from  a  publication  by 
E.  0.  Shaw  (ref  6)  and  processed  Into  a  form  usable  In  the  derived 
equations.  Shaw's  experiments  were  conducted  on  thick  barium  titanate 
disks  of  various  diameters.  His  data  were  presented  in  graphical  form 
and  consequently  some  error  In  extracting  the  data  may  be  expected. 

Table  III  extracted  from  reference  (8)  p  195,  is  Dennison 
Bancroft's  numerical  solution  to  the  complicated  classical  Pochhammer- 
Chree  solution  of  the  Isotropic  coupled  longitudinal  mode.  A.  compara¬ 
tive  table  IV  Is  obtained  from  the  longitudinal  equation  derived  In 
this  report  by  reducing  the  equation  to  the  Isotropic  form. 


i 


ft 
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Table  1.  Effect  of  length  on  Reaottant  Frequency 


Fundamental 

lgjt£th  Resonant  Frequency 

(kcV 


Sample  BT-1  (Dimeter  0.137  In.) 


0.307 

0.377 

0.348 

0.330 

0.301 

0.181 

0.161 

0.143 

0.138 

0.095 

0.063 

0.038 


363.3 

395 

333.5 

363.5 
397 
436 
491 
543 
580 
683 
763 
800 


Sample  BT--3  (Diameter  0.346  In.) 


I  0.159  390 

0.133  419 

I  435 

I  0.038  459 

I  Sample  BT"3  (Diameter  0.147  in.l 

I  0.359  226 

I  0 . 333  391 ^  3 

I  0.306  363.5 

I  0.364.  304.5 

I  337 


U*  Effect  of  Diameter  on  Resonant  Frequency  of  Tblck  Barium 
Tltanate  Disks  Reprocessed  from  Reference  6. 


(M»* 

<fW* 

(kc-ln.)* 

(VD)" 

(fD)» 

(kc«ln.) 

1.00 

5304 

0.111 

13440 

0.640 

8073 

0.095 

13710 

0.445 

10680 

0.083 

14040 

0.336 

11400 

0.071 

14840 

0.350 

11730 

0.063 

14330 

0.197 

13950 

0.055 

14530 

0.160 

13930 

0.049 

14310 

0.133 

13540 

0.044 

14100 
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Table  III.  Bancroft's  isotropic  longitudinal  cross  sectional  frequency  correc¬ 
tion  factor  as  a  function  of  Poisson's  ratio,  Mot  .  and  the  diameter  to  wave¬ 
length  ratio: 

.2 


n 


D 

2ir 

0.  1 

0.  15 

Poisson's  ratio, 

0.20  0.25 

^RL 

0.  30 

0.  35 

0.  40 

0.00 

1.  000 

1.000 

1.  000 

1.  000 

1.  000 

1.  000 

.1.  000 

0.05 

0. 99988 

0.99972 

0. 99950 

0. 99922 

0.  99888 

0. 99848 

0. 99802 

0.  10 

0. 99950 

0.99886 

0. 99798 

0. 99686 

0.  99549 

0. 99389 

0. 99206 

0.  15 

0. 99882 

0.99736 

0. 99533 

0. 99277 

0.  98967 

0.  98609 

0. 98202 

0.20 

0. 99780 

0. 99509 

0. 99138 

0. 98670 

0.  98117 

0.97480 

0. 96772 

0.  25 

0. 99632 

0.  99184 

0. 98579 

0. 97830 

0.96956 

0. 95975 

0. 94903 

0.  30 

0. 99421 

0.98728 

0. 97810 

0. 96702 

0.  95435 

0. 94049 

0. 92571 

0.  35 

0. 99114 

0. 98085 

0. 96759 

0. 95203 

0.  93486 

0.91660 

0. 89770 

0.  40 

0. 98651 

0. 97158 

0. 95310 

0. 93236 

0.  91031 

0.  88770 

0. 86503 

0.  45 

0. 97913 

0.95778 

0. 93300 

0. 90668 

0.  88003 

0.  85372 

0. 82812 

0.  50 

0. 96621 

0.  93646 

0. 90503 

0.  87383 

0.  84372 

0. 81499 

0.  78780 

Table  IV.  Cross-sectional  frequency  correction  factor  as  computed  from  the 
isotropic  form  of  the  longitudinal  frequency  equation  derived  in  this  report. 


.  ,  J^RLr  D.2_ 

D 

lU 

0,  1 

0.  15 

Poisson's  ratio, 

0.  20  0.  25 

^RL 

0.  30 

0.  35 

0.  40 

0.00 

1.  000 

1.000 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

0.05 

0. 99988 

0. 99973 

0,99950 

0. 99916 

0. 99871 

0.  99811 

0. 99733 

0.  10 

0.  99955 

0.99894 

0. 99800 

0. 99666 

O'.  99486 

0.  99246 

0. 98933 

0.  15 

0. 99889 

0.99735 

0. 99500 

0.  99167 

0. 98715 

0.  98115 

0. 97334 

0.20 

0. 99806 

0.  99576 

0. 99200 

0. 98667 

0. 97944 

0.  96984 

0. 95734 

0.25 

0. 99722 

0.99338 

0. 98750 

0. 97917 

0.  96875 

0.  95287 

0. 93335 

0.  30 

0.99603 

0.99046 

0. 98200 

0. 97001 

0. 95374 

0. 93214 

0. 90669 

0.  35 

0. 99445 

0.98676 

0. 97500 

0. 95835 

0. 93575 

0.  90575 

0. 86670 

0.40 

0. 99289 

0. 98305 

0. 96800 

0. 94668 

0. 91776 

0.  87936 

0. 82937 

0.45 

0.  99101 

0.97855 

0. 95950 

0. 93252 

0.89591 

0.87731 

0. 78406 

0.  50 

0.98890 

0.97353 

0. 95000 

0.  91667 

0. 87150 

0. 81150 

0. 73340 
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4.3  SourcejofBrror 


Basically  the  sources  of  error  are  in  the  measurement  of  di¬ 
mensions  and  frequency,  and  the  effects  of  sample-holder  pressure.  Be¬ 
cause  of  the  tolerances  of  the  ultrasonic  impact  grinder,  a  3-mll 
variation  in  diameter  occurred  along  the  length  of  some  cylinders.  As 
a  result,  this  figure  was  chosen  to  be  the  tolerance  for  length  grind¬ 
ing.  After  recording  the  frequencies  of  sample  BT-1  and  observing  the 
results,  however,  the  length  tolerance  was  reduced  to  1  mil.  Special 
attention  was  given  to  maintaining  the  holder  pressure  as  light  as  pos¬ 
sible  on  the  sample,  especially  for  long  samples,  bcess  pressure  could 
result  in  obvious  changes  in  resonant  frequency. 

No  extreme  attention  wae  given. to  the  body  temperature  of  the 
sample  other  than  making  a  rapid  frequency  measurement  at  room  tempera¬ 
ture  to  avoid  any  drift. 

Obviously,  precautions  such  as  vacuum-deposited  electrodes, 
polished  surfaces,  and  temperature  boxes  could  have  been  instituted.  It 
was  thought,  however,  that  an  accuracy  of  S  kc  was  all  that  could  be  ex¬ 
pected  from  the  instrumentation;  therefore  precise  sample  preparation 
was  not  warranted. 

B.  DiacassioN 

The  radial  and  longitudinal  frequencies  predicted  by  equations  34 
and  3S  when  plotted  against  length,  for  a  constant  diameter,  take  the 
form  shown  in  figure  2.  The  ordinate  intercept  of  the  curve  for  the 
radial  mode  is  the  ideal  frequency  fg,  and  the  abscissa  intercept  is  the 
length  value  at  which  the  quantity  CiL/D  becomes  equal  to  unity.  The 
abscissa  Intercept  of  the  curve  for  the  longitudinal  mode  is  the  length 
value  at  which  the  quantity  C2O/L  becomes  equal  to  unity. 

The  trend  of  the  theoretical  equations  in  figure  2  was  shown  to  be 
correct  by  plotting  the  measurements  made  on  sample  BT-1  in  figure  3. 

The  right-hand  portion  of  the  curve  shows  that  as  the  length  of  the  sample 
was  decreased,  the  resonance  rose  at  an  increasing  rate.  However,  the 
strength  of  the  vibration  Indicated  on  the  voltmeter  (fig.  1)  steadily  de¬ 
creased  to  a  region  of  lengths  where  two  equally  strong  resonance  readings 
were  encountered.  This  was  the  transition  region  where  the  primary  mode 
was  switching  from  the  longitudinal  to  the  radial.  Beyond  this  region,  as 
the  length  was  further  decreased,  only  single  strong  resonances  evolved, 
and  they  approached  the  ideal  fundamental  frequency  as  a  maximum  for  van¬ 
ishing  lengths. 

The  compatabillty  of  the  theory  with  the  data  in  figure  3  may  be  seen 
readily  by  noting  that  the  form  of  the  theoretical  aquations  is  linear 
when  the  squared  tents  (fD)*  and  (L/D)*  or  (OL)*  and  (D/L)*  are  used  as 
coordinate  axes,  and  that  the  intercepts  on  the  axes  are  related  to  the 
constants  of  the  material. 
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RESONANT  FREQUENCY 
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RESONANT  FREQUENCY,  KILOCYCLES/SECOND 
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Figure  3.  Changes  In  resonant  frequency  of  sample  BT-1  as  the  length 
was  ground  away.  The  curve  shows  a  similarity  to  the  theoretical 
curves  given  in  figure  2. 
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Consequently  if  the  experimental  measuremoits  are  In  agreement  with 
the  theory,  the  high  frequencies  of  sample  BT-1  plotted  on  the  axes  (reso¬ 
nant  frequency  x  diameter)^  and  (L/D)^  should  form  a  straight  line.  The 
same  argument  applies  to  the  low  frequencies  of  sample  BT-1;  when  these 
data  are  plotted  on  the  axes  (resonant  frequency  x  length)^  and  (D/L)^,  a 
straight  line  should  also  result  Figures  4  and  5  show  that  this  is  in¬ 
deed  the  case..  Assuming  then  that  the  theory  Is  correct,  the  equations 
may  be  fitted  to  the  experimental  data  and  the  values  of  the  constants 
may  be  calculated  from  the  intercepts 

G  =  111.3  In.-kc  H  =  82.2  In.-kc 

0.7C7  C^=  0.370 

Comparative  values  of  G  and  H  are  available  In  reference  1. 

G  =  111.7  in.-kc  H  =  82.7  In.-kc 

No  literature  Is  available  to  provide  direct  comparison  with  the 
and  C^  constants.  However,  the  components  of  the  constants  are  readily 
available: 

\i  =  0.3  (ref  7) 

R  =  2.049  (refU) 
o 

Inserting  these  values  into  equations  32  and  33  yields 

Cj  =  0.726 

C_  =  0.356 

which  are  in  good  agreement  with  the  values  0.757  and  0.37  obtained  ex¬ 
perimentally  . 

As  further  evidence  that  the  constants  Cj^  and  C2  are  Independent  of 
dimensions,  the  measurements  obtained  from  samples  BT-2  and  BT-3  were 
plotted  In  figures  6  and  7  together  with  the  fitted  equations  obtained 
from  figures  4  and  5;  the  degree  of  the  fit  Is  sufficient  to  warrant  the 
conclusion  that  the  constants  are  Independent  of  dimensions. 

In  all  experimental  work  reported  herein,  the  cylinder  length  was 
the  dimension  varied.  In  an  effort  to  verify  the  theory  for  diameter 
changes  as  well,  the  data  presented  by  E.  0.  Shaw  In  reference  6  were 
examined  for  conformance  to  the  radial  theory.  These  data  were  published 
In  graphical  form  so  some  error  In  extracting  the  data  was  expected;  how¬ 
ever,  when  they  were  plotted  in  figure  8,  they  followed  the  form  of  the 
theoretical  radial  equation.  The  intercepts  differed  from  those  obtained 
experimentally,  and  the  conclusion  is  that  a  different  barium  titanate 
formulation  was  used  by  Shaw. 
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RESONANT  FREQUENCY  -  DIAMETER  PRODUCT  SQUARED 
(KILOCYCLES/SECOND  X  INCH)^ 


(L/D)^ 


l''ijiurc  ‘1.  IlirJ'  I’l'eqiu'ncy  data  ohtaiiu'd  I'l-oni  samjjle  B  I'-l  illiistralt  s 
till!  taiidi'iicy  of  this  data  to  follow  llu!  linear  requirements  of  the  radial 
equation  .  Fitting  the  radial  eqtmtion  to  the  data  forces  the  ordinate 
interceril  to  yield  the  constant  and  the  abscissa  intercept  the  constant 

OIg/-. 


RESONANT  FREQUENCY  -  LENGTH  PRODUCT  SQUARED 
{KILOCYCLES/SECOND  X  INCH)* 


5000  t  1  — 1  - ^ - i  i  i  1  I  I  > 
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(D/D* 

Figure  5.  I.ow  frequency  data  obtained  from  sample  BT-1  illustrates 
the  tendency  of  these  data  to  jneet  the  linear  requirements  of  the 
longitudinal  theory.  l'’ittiii';  the  longitudinal  equation  to  tlie  data  forces 
the  ordinate  intercept  to  represent  the  constant  and  the  abscissa 
intercciH  the  constant  (l/C2)'^. 
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RESONANT  FREQUENCY-  DIAMETER  PRODUCT  SQUARED 
{KILOCYCLES/SECOND  X  INCH)^ 


(L/D)^ 


Figure  6,  Plot  of  data  for  the  radial  mode  specimens  prepared  from 
sample  BT-2  (identical  material  to  BT-1  but  of  a  larger  diameter)  and 
the  curve  determined  from  BT-1  specimens  (Figure  4).  The  compar¬ 
ison  shows  that  the  constants  G  and  C  are  representative  of  the  material 
and  independent  of  cylinder  dimensions. 
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RESONANT  FREQUENCY-  LENGTH  PRODUCT  SQUARED 
(KILOCYCLES/SECOND  X  INCH)^ 
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Figure  7.  Plot  of  data  for  the  longitudinal  mode  specimens  pi’epared  from 
sample  BT-3  (identical  material  to  BT-1,  slightly  larger  diameter)  and 
the  curve  determined  from  BT-1  specimens  (figure  5).  The  comparison 
shows  that  the  constants  H  and  C2  are  representative  of  the  material  and 
independent  of  the  cylinder  dimensions. 


RESONANT  TREQUENCY-  DIAMETER  PRODUCT  SQUARED 
(KILOCYCLES/SECOND  X  INCH)* 


LENGTH-DIAMETER  RATIO  SQUARED 
(L/D)* 


Figure  8.  Reworked  d^ta  from  ref  6  for  multidiameter  cylinders  also 
shows  excellent  conformity  with  the  linear  requirements  of  the  radial 
equation. 


26 


Limited  data  were  presented  in  reference  8,  but  they  allowed  a  test 
of  the  modified  Rayleigh  theory  on  isotropic  steel  cylinders.  Using  the 
material  properties  given  in  the  report,  the  constants  Cg  and  H  were  cal¬ 
culated,  and  the  resulting  longitudinal  equation  accurately  predicted  the 
resonant  frequencies  obtained  experlmentallv  . 


However,  in  all  fairness  it  should  be  mentioned  that  more  extensive 
data  on  the  longitudinal  vibration  of  steel  rods  was  published  in  ref¬ 
erence  9,  page  148,  which  did  not  agree  completely  with  the  theory.  In 
that  work  the  claimed  values  of  material  properties  did  not  produce  cor¬ 
rect  values  of  H  and  C2,  although  the  linear  requirements  were  met.  The 
possibility  that  the  steel  was  slightly  anisotropic  was  recognized  in 
the  report,  and  is  accepted  here  as  the  reason  for  deviation  from  the 
isotropic  theory.  Additional  information  was  unavailable  but  necessary 
in  order  to  employ  the  anisotropic  equation. 


Dennison  Bancroft's  numerical  solution  of  the  classical  isotropic 
Pochhammer-Chree  equation  of  the  coupled  longitudinal  mode  is  listed  in 
table  III.  According  to  this  data  the  linear  requirements  of  the  iso¬ 
tropic  longitudinal  equation  are  contradicted.  Table  IV  shows  that  only 
a  slight  error  exists  in  frequency  prediction  between  the  two  methods, 
but  this  error  will  naturally  be  Increased  as  the  dimension  ratio  approaches 
unity. 

Figure  9  la  a  brief  comparison  of  the  two  tables  and  the  discrepancy 
cannot  be  explained.  It  will  only  be  noted  that  in  the  region  of  p  s  0.38 
the  two  tables  are  in  closest  agreement,  that  this  is  the  region  of  both 
the  steel  and  piezoelectric  ceramic  data,  that  the  equations  developed  in 
this  paper  have  the  advantage  of  being  able  to  handle  some  anisotropic 
materials  vibrating  in  either  mode,  and  that  all  the  experimental  data  pre¬ 
sented  in  this  report  do  appear  to  follow  the  linear  requirements  of  the 
theory • 


6.  CX)WCLUSION 

The  experiments  performed  to  verify  the  equations 
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showed  that  the  trend  of  the  equations  is  correct,  the  constants  are  in 
agreement  with  those  in  the  literature  for  barium  tltanate  and  steel,  the 
equations  are  unrestricted  by  anisotropic  effects  of  piezoelectrics,  and 
the  equations  are  valid  for  any  dimensions. 
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LONGITUDINAL  CORRECTION  FACTOR 


Figure  9.  Comparison  of  Bancroft's  iminorical  solution  of  coupled  longitudinal 
vibrations  with  the  i;;<iirt)pic  form  of  the  longitudinal  equation  derived  in  this 
report. 
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Theso  equations  fill  a  void  that  has  hindered  the  design  of  cylin¬ 
drical  microminiature  ceramic  transformers  and  filters.  The  derivations 
themselves  avoided  the  use  of  higher  mathematics  and  tedious  tensor 
equations.  Their  simplicity  opens  the  door  to  the  derivation  of  addition¬ 
al  design  equations  that  eliminate  the  need  for  the  time-consuming  process 
of  experimentallydetermlning  an  empirical  relation. 

In  the  field  of  piezoelectrics,  the  most  valuable  form  of  the  fre¬ 
quency  equations  is: 

cr  *  ' 

Cf  T  *  » * 

The  constants  C^,  G,  and  H  should  subsequently  be  listed  in  the  lit¬ 
erature  for  various  materials,  strengths  of  polarization,  and  tempera¬ 
ture. 
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SYMBOL  GLOSSARY 


D 

E. 


In 

W 

P 

AW 

a 


a 


6 

6 

6 

P 


L 

* 

R 


Equation  No. 


Grouping  of  radial  material  constants,  dimensionless  32,  19 

Grouping  of  longitudinal  material  constants,  dimensionless  33,  23 

Diameter,  in.  1 

Longitudinal  modulus  of  elasticity,  Ib/in.*  2 

Radial  modulus  of  elasticity,  Ib/in.^  1 

Actual  radial  resonant  frequency,  kc  34 

Ideal  radial  resonant  frequency  for  disk  with  zero  thick¬ 
ness,  kc  1 

Net  longitudinal  force,  lb  6 

Net  phantom  force,  lb  12 

Net  radial  force,  lb  9 

Radial  frequency  constant,  in.-kc  10 

Longitudinal  frequency  constant,  kc  7 

Longitudinal  spring  constant,  Ib/ln.  6 

Radial  spring  constant,  Ib/ln.  9 

Maximum  kinetic  energy,  in. -lb  15 

•Length,  in.  2 

Effective  lumped  mass  of  cylinder  in  longitudinal  vibration, 
lb/sec®/ln.  7 

Effective  lumped  mass  of  cylinder  in  radial  vibration, 

Ib/sec^/ln.  10 

Longitudinal  constant  of  the  phantom  force,  dimensionless  22 

Radial  constant  of  the  phantom  force,  dimensionless  12,  30 

Input  work,  in. -lb  11 

Phantom  work,  in. -lb  13 

Work  available  lor  radial  vibrations,  inrlb  14 

Actual  longitudinal  resonant  frequency,  kc  35 

Ideal  longitudinal  resonant  frequency  for  rod  with  zero 
diameter,  kc  2 

Maximum  longitudinal  deflection,  in.  5 

Intermediate  radial  deflection,  in.  15 

Maxlmvun  radial  deflection,  in.  8 

Density,  lb-sec*/in.*  7,  10 
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SYMBOL  GLOSSARY  (Cont'd) 


Equation  No. 


^L  Longitudinal  strain,  dimensionless  26 

J,  J  Bessel  functions,  ref.  S,  p  494  1 

1,  o  *  ’ 

Poisson's  ratio  associating  radial  and  tangential  strains  in 
the  radial  plane,  dlMnslonless  1 

11  Isotropic  Poisson's  ratio,  dimensionless  3,  4 

Poisson's  ratio  relating  the  radial  strains  to  the  longitudinal 
strains,  dimensionless  26 

Function  defined  In  equation  (1),  dimensionless  1 

t  Time,  sec  15 

Y_  Numerical  constant,  dimensionless  28 

R  ' 
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